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ABSTRACT 


The goal of this project was to develop computational tools for the calculation of 
electronic structure of molecules containing heavy atoms, and to use these tools in the 
study of catalytic processes, with the overall objective of gaining an understanding of 
the catalytic process which could be used to design more efficient catalysts. The main 
catalytic system of interest was the combustion of hydrogen on platinum surfaces. Un- 
der this project, a flexible Dirac-Hartree-Fock (DHF) program has been developed, a 
code to calculate DHF correlation energies at the second-order Mpller-Plesset perturba- 
tion (MP2) level is almost complete, ]and code to include correlation at the MCSCF and 
MCSCF/MP2 level is planned. The tools so far developed have been validated and used 
to calibrate some more approximate methods, and applied to investigate the importance 
of relativistic effects in the bonding of hydrogen to platinum. 


The chemistry of the heavy elements is strongly affected by relativity, and any ap- 
proach to computing properties must therefore be based on a relativistic theory. Because 
of the increased difficulty of relativistic molecular structure calculations over nonrela- 
tivistic calculations — one to two orders of magnitude in storage, for example ~ it has 
been usual in the past to make approximations. Two common approximations are rel- 
ativistic effective core potentials (RECPs), in which all the relativistic effects are incor- 
porated into a potential which also replaces the core electrons, so that only the valence 
electrons are treated explicitly by the theory, and first-order perturbation theory (PT) 
using only the mass- velocity and Darwin terms (MVD) of the one-electron first-order 
Hamiltonian, which retains all the electrons, but truncates an expansion of the full rel- 
ativistic Hamiltonian so that the approximate perturbed Hamiltonian can no longer be 
used variation ally. However, neither of these methods is fully satisfactory, and it is nec- 
essary to ensure that they are accurately reproducing the relativistic effects. Thus efforts 
have been made in this project to develop methods based on the Dirac equation which 
will enable accurate relativistic molecular calculations to be made and calibrate the more 
approximate methods. 

The first stage was to develop a Dirac-Hartree-Fock (DHF) code for polyatomic 
molecules. Because the field is a new one, there were many issues to be considered 
and explored in the implementation. These issues are discussed at length in two pub- 
lications [1,2]. Several criteria were used to shape the code development. As a first 
generation code, it was desirable to get a code running as far as possible by adaption 
of existing codes or techniques. Thus, the integrals were taken from an adaption of 
MOLECULE [3]. Second, it was assumed that disk space limitations would preclude the 
sorting or transformation of integrals. Consequently, the program was designed to con- 
struct Fock matrices from the unordered scalar integrals which are directly output from 
MOLECULE. The use of MOLECULE also dictated the implementation of symmetry 
via symmetry-adapted basis functions rather than the petite list approach [4]. One im- 
portant feature of the code is that, despite the use of scalar basis functions, the internal 
representation of the basis is a 2-spinor basis. This averts many of the problems of linear 
dependence of the basis in the only other code in existence at the time of development, 
namely MOLFDIR [5]. 

After the initial development of a closed-shell DHF code, other features were 
added to the program to increase its flexibility. These included convergence accelera- 
tion using the DIIS algorithm [6], the evaluation of molecular dipole and quadrupole mo- 
ments, the application of a finite (constant) electric field, and various spinor population 
analyses. Subsequently, the code has been extended to allow calculations on molecules 
with up to two electrons in open-shells, including a reference function consisting of two 
closed-shell configurations. 


Some of the issues surrounding the use of the DHF method were investigated, 
such as the question of whether restricted or unrestricted kinetic balance was to be 
preferred [7], and the appearance of intruder states which appeared in some of the cal- 
culations [8]. The second of these was prompted by the investigation of the neglect of 
some integrals whose contribution was argued to be small, done in an effort to reduce 
the space demands of DHF calculations. 

Having developed a molecular code it was necessary to have some means of gener- 
ating basis sets. An adaption was made to the program GRASP [9] to enable generally 
contracted atomic basis sets to be developed and printed in a format suitable for input 
to MOLECULE. The adaption has been extended to permit exponent optimization [10] 
using gradient techniques, a feature which may be important for heavy elements. 

Several sets of calculations have been done with the DHF program. In a series 
on the group IV elements, the properties of the tetrahydrides, dihydrides and monoxides 
have been calculated, and comparisons made with RECP and PT calculations [11-13], 
These calculations revealed deficiencies in some of the published RECPs, and suggest 
that a broader calibration study would be profitable. Two other studies have produced 
estimates of relativistic effects in alkali metal compounds [14,15]. 

In a first step towards catalytic studies, the DHF method has been applied to 
the mono- and dihydrides of platinum [16]. Comparisons were made with PT and the 
use of the second-order Douglas-Kroll transformed Hamiltonian developed by Hess [17]. 
The last of these should give a fairly accurate picture of the spin-free relativistic effects 
in the valence shell, and enable an estimate of the importance of spin-orbit effects. (It 
is primarily because of the spin-orbit interaction that relativistic calculations are much 
more expensive than nonrelativistic calculations.) The calculations showed that even for 
a closed shell molecule the spin-orbit effect could have a significant effect on molecular 
properties. 

The second-order Douglas-Kroll transformed Hamiltonian is an approximation 
to the Dirac Hamiltonian which attempts to decouple the large and small components 
of the wave function, and then separate the spin-orbit and spin-free terms. It is in the 
first part of the process that the approximations are made. If the small component is 
retained, it is still possible to separate the spin-orbit and spin-free terms, in a rigorous 
manner. This has been achieved and exploratory calculations on atoms performed which 
demonstrate the potential of the use of the spin-free terms of the Dirac equation [18]. 
The implications of this separation are far-reaching. The cost of a correlated calculation 
with the spin-free Dirac Hamiltonian is little more than that of a nonrelativistic calcula- 
tion. It would thus be possible to use this method in conjunction with the DHF method 
at a lower level of calculation to determine whether spin-orbit effects were important, 
and if not, the higher level calculations may be done much more cheaply with the spin- 


free Hamiltonian. 

Another major area of program development which is in progress is the transfor- 
mation of integrals from the scalar basis in which they are generated to the 2-spinor ba- 
sis employed in parts of the DHF program, and thence to supermatrix form. The reason 
for these developments is that the DHF program, while written efficiently for the circum- 
stances under which it was expected to be used, is not the most efficient implementation 
possible, and with the rapid increase in available computer resources over the years, it is 
now possible to consider transforming the integrals. The use of ordered integrals, either 
in the scalar basis or in the 2-spinor basis, would considerably speed up the construction 
of the Fock matrix, and even more so if supermatrices were constructed. Furthermore, 
in order to proceed beyond the SCF level and include electron correlation it is neces- 
sary to transform the integrals to the molecular 4-spinor basis. Therefore, a considerable 
amount of effort has been spent on analyzing the integral ordering and transformation 
for the DHF problem. Much of this work was used in preparation for the NATO Ad- 
vanced Summer Institute in Vancouver, BC, in August 1992 [1], 

The transformation of the integrals is part of the second stage of the code devel- 
opment, which is to include methods for electron correlation, a necessity for accurate 
description — and in many cases even for qualitative description - of metal compounds. 
The first step in the inclusion of electron correlation has been made at the simplest level, 
namely second-order M0ller-Plesset (MP2) perturbation theory, which has been devel- 
oped for systems involving up to two open shells [19], following recent open-shell nonrel- 
ativistic theories [20]. The method has been applied to the coinage metal hydrides [21]. 
The coding of the sections for two open shells is nearing completion. 

M0ller-Plesset (MP2) perturbation theory is useful where the wave function is 
dominated by a single configuration. This is much less frequently the case in the (rela- 
tivistic) double-group representation. Therefore it is essential to develop multi-reference 
correlation methods. To this end, a collaboration has been initiated to develop a MC- 
SCF code based on the SIRIUS code [22]. At this point, most of the theoretical ground- 
work has been completed [23]. The use of MP2 theory based on a MCSCF reference has 
been shown to be useful [24]; thus the development of an MP2 code was a preliminary to 
the development of an MCSCF MP2 code. 

In summary, a flexible DHF program has been developed, a code to calculate 
DHF correlation energies at the MP2 level is almost complete, and code to include cor- 
relation at the MCSCF and MCSCF/MP2 level is planned. The tools so far developed 
have been validated and used to calibrate some more approximate methods, and applied 
to investigate the importance of relativistic effects in the bonding of hydrogen to plat- 


mum. 
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